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ABSTRACT 



Context. With the discovery of high-energy y-ray emission from the Andromeda galaxy (M31) by the FerinifLAT collaboration, 
normal galaxies begin to arise from the shadows for the first time, providing insight into cosmic ray acceleration in external galaxies. 
Aims. We search for high-energy y-ray emission from those galaxies in the Local Group that have so far not been investigated: M 81, 
M 83, IC 342, Maffei 1, Maffei 2, and M 94. 

Methods. FennifLAT public data from August 4, 2008 to January 1, 201 1 were analysed for these galaxies. We compared the results 
to other starburst and normal galaxies detected so far at high energies: the Magellanic clouds, M 31, and the starburst galaxies M 82 
and NGC 253. 

Results. No significant detection is found in the data for the sources in our sample, and we derive upper limits on their photon flux. 
After comparing the results to other Local Group objects, we find that the derived upper limits are fully compatible with expectations 
from cosmic rays interacting with the interstellar medium within the host galaxies. In the case of M 33 and M 83, a detection in 
FermifLKI data should be imminent. The expected fluxes for the other sources in the sample are below the sensitivity of Fermi/LAT, 
even after 10 years of observation. Collective emission from compact objects in the host galaxies is also found to be negligible 
compared to the expected emission from cosmic ray interactions. 

Key words, gamma rays: galaxies ~ ISM: cosmic rays ~ radiation mechanisms: non-thermal 



1. Introduction 

The usual suspects for extragalactic high-energy y-ray emis- 
sion are active galactic nuclei (AGNs), and more particularly 
sources presenting powerful, relativistic jets pointing close to 
the line of sight, the so-called blazars, which are BLLac ob- 
jects or flat spectrum radio quasars. Their contribution to the 
extragalactic diffuse y-ray emission as unresolved sources has 
been extensively discussed (see e.g. IStecker & Salamonlll996t 
Dermer 2007). However, less powerful sources begin to be re- 
vealed as high-energy emitters. For instance, the starburst galax- 
ies M 82 an d NGC 253 have be en detected at high energies with 
FermifLNV jAbdo et al.ll2010ah . Recently, we have proposed the 
association of a previously unidentified Fenni/LAI source from 
the FermifLPS first source catalogue (IFGL, Abdo et al. 2010b) 
with the archetypal Seyfert2 galaxy NGC 1068 ( Lenain et al] 
12010 ). We also discussed there the origin of the high-energy y- 
ray emission of NGC 4945, another Seyfert 2 galaxy encompass- 
ing a starburst region in its core, and argued that in this case, the 
starburst activity within may be the dominant emitter with re- 
spect to the central AGN, while the AGN activity might play the 
main role at high energies in the case of NGC 1068. 

Normal galaxies also begin to arise from shadows in the GeV 
range, as can be seen by the recent discovery of y-ray emission 
with Fermi/LAT from our neighbour, the Andromeda galaxy 
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rMSl. lAbdoeTal ] |2010d) . Normal and Seyfert galaxies have 
been thought to contribute to a s mall amount o f the extragalac- 
tic diffuse y-ray emission (see e.g. lStrong et al.i ri976). However, 
recent studies have shown that their contribution can in fact be 
domi nant compared to the one from blazars (see e.g. Fiel ds et aP 
I2010t) . We present here a study of the search for high-energy 
y-ray emission from the major galaxies pertaining to the Local 
Group. 

We first describe the sample of sources and the method for 
the high-energy y-ray data analysis in Section |2l before dis- 
cussing our results in Section [3] and drawing some conclusions 
in Section|4] 



2. Observations and data analysis 

We focus our study on the major galaxies of the Local Group 
(see e.g. Karachentsev 2005) for which no results at high en- 
ergies have been discussed so far, that is M81, M83, IC342, 
Maffei 1, Maffei 2, and M 94. The other galaxies from the Local 
Group for which high-energy em ission has been repo rted are the 
smafl Magellanic cloud (SMC, Abdo et al.' 2010d), the lai-ge 
Mageflanic cloud (LMC, Abdo et al. 2010e), the Andromeda 
galaxy (M31, Abdo et al. 201 0^), M 82 (lAbdo et al.l l201()ah . 
Cen A (NGC5128.iAbdo et aljl2010jigh. and NGC 253 l Abdo 
et al. l2010al) . Our primary goal is to search for signatures of 
high-energy y-ray emission from interactions of cosmic rays ac- 
celerated by stellar processes (e.g. supernovae) with the inter- 
stellar medium in the host galaxy, as opposed to y-ray emission 
from AGN activity. 
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For the different sources in our sample, the FermifLAI data 
analysis was performed using the publicly available data from 
August 4, 200 8 to January 1, 20 11. We used the unbinned likeli- 
hood method (lAtwood et al.i 2009) from the public Science Tools 
(version v9r 18p6) analysis software provided by the Fermi col- 
laboration. The data analysis was carried out following the pro- 
cedure recommended by the Fermi/LAH^ collaboratiorQ, using 
diffuse class events, with the P6_V3 instrument responses, in the 
energy range of 200MeV-200GeV. The Galactic diffuse emis- 
sion was modelled using the public model gll_iem_v02, while 
the isotropic model isotropic_iem_v(52 was used to account 
for both extragalactic, unresolved diffuse emission and residual 
instrumental background. These two models are available as part 
of the public Science Tools. For each source in the sample, we 
selected events from a circular region of 10° of radius around the 
nominal position of the galaxy of interest. Given the distances of 
the sources considered for our analysis - the closest one being 
Maffei 2 lying at 2.8 Mpc - the objects are assumed to be point- 
like at high energy. Indeed, the apparent dimensions of the host 
galaxies are mostly below 10', i.e. smaller than the point spread 
function of the Fenni/LAT in the major part of the considered 
energy range. 

In the modelled reconstruction of the sources, all t he ob- 
jects included in the IFGL catalogue (lAbdo et al iHMob) are 
accounted f or. Moreover, other sources, for which test statistics 
(TS, see e.g.' Mattox et al.iri996h are above 25 - roughly equiva- 
lent to a 5cr detection - and which are not reported in the IFGL 
catalogue, are also included in the source models. These poten- 
tial source candidates were selected from the Ferm//LAT count 
maps above 1 GeV, where the Ferm//LAT point spread function 
is below ~0.5°. For each object added in this way to the model, 
as well as for the source of interest, we first run a likelihood test, 
using gtlike, to assess whether the new introduced source is 
significant. If the corresponding TS is above 25, revealing a de- 
tection, the position of the new source is optimised using the tool 
gtfindsrc, whereas for the sources included in the IFGL cat- 
alogue, we fixed their position to the values from the catalogue. 
The spectral parameters of the additional source are then refined 
on the corresponding best-fit position, using gtlike. 

For all the sources of interest in the sample, the correspond- 
ing TS are found to be below 25. We thus computed upper limits 
at 2cr confidence level on their flux in the 200 Me V-200 Ge V 
energy range. No specific energy spectra are assumed to derive 
these upper limits, the photon indices of the sources of interest 
are left free to vary in the models. Table [T] summarises the re- 
sults of our analysis on the source sample, as well as results for 
other galaxies studied in the Fermi era. Following are some notes 
concerning peculiar sources. 



M81 



For the search of y-ray emission from M 81, the strong starburst 
emitter M82 is located 0.62° away from M81, hence contami- 
nating the region around M 81, given that the point spread func- 
tion of the LAT degrades strongly at low energies, which pre- 
vents us from deriving constraining upper limits on the flux of 
M81. The corresponding upper limit given in Table [T] should 
thus be considered as conservative. 




Fig. 1. Smoothed count map of the region of interest around 
M83, in the IGe V-200 GeV energy range. The nominal po- 
sitions of M83 and the blazar 2E3100 are shown in white. 
The best-fit position, at 68% and 95% confidence levels, of the 
Fermij'LPS excess is shown in blue. 



M83 

A significant signal is found in the vicinity of M83, with a 
TS=30.1 (5.5cr). However, the best-fit position of the excess 
turns out to be compatible with the nominal position of a 
blazar, 2E3100, lying at 0.33° from the position of M83. The 
FermifLAI source is not spatially compatible with M 83, as can 
be seen in the count map shown in Fig.[T]for the 1 Ge V-200 GeV 
energy range. 

IC342 

IC 342 lies at 10.6° from the Galactic plane. Given the proximity 
of the Galactic plane for this source lying at b = -1-10.6 deg, we 
extracted the events from a region of 10° of radius centred on 
a position displaced by 5.7° towards north-east compared to the 
nominal position of IC 342, in order to reduce the fraction of the 
emission from the Galactic background in the region of interest. 
For this analysis, four sources were added to the source model 
in addition to those from the IFGL catalogue, and the signal 
extracted around IC 342 amounts to 1.3cr. 



Maffei group 

The two galaxies Maffei 1 and Maffei 2 are separated by only 
0.68°, making it difficult to derive high-energy properties for 
these two objects separately with Fermi/LAJi, given the angular 
resolution of the LAT0. Moreover, they are located at less than 1° 
from the Galactic plane. For this reason, no precise constraints 
on their supernova rate or the total gas mass were found in the 
literature. 



M94 

M94 is a spiral galaxy viewed almost face-on. The closest 
source from the IFGL catalogue lies at ~4° from M94, eas- 



' see |http://fermi.gsfc.nasa.gov/ssc/data/analysis| 



see |http://www-glast.slac.stanford.edu/software/IS/glast-lat-performance.hti 
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Table 1. Properties of galaxies from the Local Group. 



Source 


d 








TS 




kpc 




lO'M© 


lO'*^ ph s"' 




Milky Way 




0.02 ±0.01*" 


6.5 ± 2.0<-' 


0.12 ± 0.03*3* 




LMC 


SO*"' 


0.005 ± 0.002<^> 


0.67 ± 0.08'*' 


(7.8 ± 0.8) X 10^3'^' 


1110.7™ 


SMC 


61<*' 


(1 +0.2) X lO"'*" 


0.70 ±0.06'"" 


(1.6 ± 0.4) X 10^3(11) 


138.8"" 


M31 


780(12) 


0.045 ± 0.015 


7.66 ±2.38"^' 


(6.6 ± 1.4) X 10"^<"" 


28.9<'*' 


M33 


847(15) 


0.050 ±0.015 


2.23 ±0.84'"'' 


< 5.0 X 10"^"*' 


< 25'"** 


Matfei 2 


2800'"' 


7 


7 


< 1.66 


0.0 


Maffei 1 


3010*"*' 


7 


7 


< 3.03 


2.5 




'!non( 19) 
5Zo\r ' 


U. to ± U. iU ' 


/I n _i_ n o(21) 
4.U ± U.o 


< D.4U 


1 O 


M82 


3530<22) 


0.2±0.1<'> 


2.5±0.7<" 


2.39±0.75<'> 


46.2<" 


NGC4945 


3600<23' 


0.3 ± 0.2<^^' 


4 2(25) 


4.09 ± 0.92<2« 


85.3(2«) 


M81 


3630<"> 


0.008 ± 0.002<28' 


5.16± 1.72<2'') 


< 1.55 


0.0 


NGC253 


3940(30) 


0.2 ±0.1*" 


2.5 ± 0.6<" 


1.12±0.78<'> 


23.0<" 


M83 


4500*"' 


0.050 ± 0.025*'" 


5.5 ±1.1(3" 


< 1.33 


5.8 


M94 


4660<32) 


0.04 ± 0.02<"> 


0.56 ± 0.1 1<"> 


< 3.19 


0.0 


NGC 1068 


WOOO*''" 


0.20 ± 0.08"^> 


4 4(36) 


33.8 ± 8.9<2«' 


68.6P«) 



References. (l)lAbdo et al' |l2010al); (2) 'Damd (11993'); (3)IStron£ et al.l (120101); (4) IPietrzvnski et al.l (12009'); (5)'Tammann et al.' (11994)); (6) 
IStave lev-Smithetal. (2003); Fukui et al. (2008); (7) Abdo et al. (2010e); (8) Hilditch et al. (2005); (9) Seward & Mitchell (1981); (10) Stanimirovic 
et al. (il999l) ; iLerov et al.i ( 1 20070; (IDiAbdo et aL (i2010d); (12) Stanek & Gamavich ( 1998); (13) Nieten et al. (2006); Braun et al. (2009); (14) 
Abdo et al.i (l2010c(); (15)lGalleti et all (l2004ft; (16 ) Gratier et al. (2010); (17) Karachentsev et al. (2003a); (18) Finaerhut et al. (2003); (19) 
Karachentsev et al. (2004); (20) Ando et al. (2005); (21) Rots (1979); (22) Karachentsev et al. (2002a); (23) Karachentsev et al. (2002b); (24) 
Lenc & Tinaav (2009); (25) WeiB et al. (2008); (26) Lenain et al. (2010); (27) Karachentsev & Kaisin (2007); (28) Matonick & Fesen ( 199j); 
(29)jGurovich et al. (2010); (30) Karachentsev et al. (2003b); (31) Lundg ren et al.l (12 004); (32) Karachentsev et al. (2003c); (33) Bagetak os et^ 
(1201 Ih ; (34) ,Gil de Paz et aH2007.) ; (35) .Wilson & Ulvestad C1982.) ; .Blietz et al.l (fl994' ); .Mannucci et al., (.2003,) ; (36) Sage et al., a99a) . 



ing any analysis of this region. The signal extraction from M 94 
amounts to O.Ocr. 



3. Discussion 

High-energy y-ray emission is expected from starburst galaxies, 
as well as normal galaxies, through cosmic ray interactions with 
the interstellar medium. Cosmic rays, which are mainly made up 
of charged, light hadrons and leptons that are accelerated in these 
galaxies, interact with the ambient interstellar medium, creating 
pions which will in turn decay into y-rays for the hadrons, and 
bremsstrahlu ng and inverse Compton radiations for the leptons 
(see e.g. Stro ng et al.ll2010l). Based on a simple model. Pavlidou 
& Fields (200l|) suggest there is a relationship between the ex- 
pected y-ray luminosity, the supernova rate 7?sn. and the total 
gas mass Mgas for a given galaxy; 



Fy(> 100 Me V) =2.34 X 10" 



r>MW 



M, 



gas 



o 



100 kpc 



(1) 



ph cm 



where is the supernova rate in the Milky Way and d the 
distance of the considered galaxy. 

Figure |2] shows the y-ray luminosity measured with 
FermifLP^ for the different known high-energy emitting star- 
burst and normal galaxies, against /?sn xMgas, and the upper lim- 
its at 2cr confidence level on the fluxes derived from the sources 
in our sample. For comparison, we added the Seyfert 2 galax y 
NGC 1068, also detected at high energies dLenain et alJl2010l) . 
This object enco mpasses a starburst re gion in its centre as well, 
but we argued in iLenain et alJ (1201 Oh that the associated high- 
energy y-ray emission is more likely due to the central AGN ac- 
tivity than from the starburst activity only, owing to its relatively 
high y-ray luminosity. 



,45 
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'w 



-10" 



10" 



10" 



10" 



10" 



__\__in' 



SMC 

t 



LMU 



10" 



10^ 



10'^ 

RsN X 



10^ 



10^ 
yri) 



Fig. 2. y-ray luminosity, or the upper limit at 2cr confidence level 
on the luminosity, of the different sources in the sample and 
other known high-energy emitting starburst and normal galax- 
ies, given against the supernova rate times the total gas mass in 
these objects. The red square points show the expectations on 
the luminosity from the model of Pavlidou & F ields. (2001.) . ac- 
counting for the uncertainties on /?sn and Mgas- 



As can be seen in Fig. |2] the upper limit reported by the 
ferm;/LAT collaboration on M33 (Abdo et al. 2010c) is very 
stringent compared to expectation, if the high-energy flux is in- 
deed directly connected to R^n x A/gas, that is, if the high energy 
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emission is solely due to cosmic ray interactions with the in- 
terstellar medium. In such a case, one could very soon expect 
a detection of this source with the FermifLPX. A non-detection 
would mean a different cosmic ray content in this galaxy, e.g. 
a lower supernovae rate than previously thought or a particu- 
larly efficient electron escape from the host galaxy to prevent 
high-energy emission. It is also noticeable that the current upper 
limit on M83 reported here also shows that a positive detection 
could be imminent. Apart from the case of NGC 1068, where the 
high-ener gy y-ray emis sion could be dominated by the central 
AGN (see lLenainet"al]|2010h . the observatio ns are fully consis- 
tent w ith the expectations from the model of iPavIidou & Fieldsl 
(I2001h . 

An alternative expectation of y-ray emission from these 
sources resides in the recent discovery by iDobler et al.l (1201 Ol) 
of giant y-ra y bubbles in th e inner Milky Way, whose origin is 
still unclear Su et al.l (1201 Ot) extensively discuss different possi- 
ble mechanisms for the origin of this large-scale y-ray emission. 
This emission more likely comes from inverse Compton scat- 
tering of a hard population of cosmic ray electrons, and may 
have the same origin as the microwave haze seen in WMAP data 
dPobler & Finkbeiner 2008). If such high -energy emission is 
common in normal galaxies, one could also expect to detect that 
signal in the galaxies of the Local Group, in addition to the emis- 
sion discussed above. In this case, the high-energy y-ray flux of 
these objects i s expected to be la r ger th an the flux predicted in 
the model by iPavlidou & Fieldsl (l200lh . In the particular case 
of M33, such emission can already be ruled out given the con- 
straining upper limit put on its y-ray flux. 

High-energy y-ray emission could also arise from the collec- 
tive emission of galactic sources within the host galaxies, such 
as pulsars or pulsar wind nebulae. To evaluate the contribution 
of such emission, we used the luminosities of pulsars and the 
few pulsar wind nebulae detecte d by ferm//LAT (Abd o et al.l 
1201 Oht lAckermann et al.l 1201 ll) . in order to estimate the ex- 
pected high-energy flux from the collective emission, assuming 
that these Galactic objects are located in a hypothetical galaxy 
located at 3 Mpc, which is about the mean distance for the source 
in our sample. We intentionally over-estimate the predicted flux 
by assuming 1000 of such sources in the host galaxy, which 
is much more numerous than what is currently detected in our 
own Galax)|3. Even under these highly optimistic conditions, 
the overall contribution from pulsars is expected to be F{E > 
100 Me V) ~ 5 X 10"''* erg cm"- s ', and the one for the pulsar 
wind nebulae amounts to F{E > 100 MeV) ~ 10"'"' erg cm"^ s"' . 
Such estimations are well below the sensitivity performances of 
the FermifLPS within two years of observation^ The total con- 
tribution of galactic objects in the sources of our sample is thus 
expected to be negligible compared to the global emission ex- 
pected from cosmic ray interactions with the local interstellar 
medium. 



4. Conclusions 

We have searched for high-energy y-ray emission from galax- 
ies belonging to the Local Group. Apart from the Andromeda 
galaxy (M31), already established as a source of y-ray emis- 
sion by the Ferm ifLAT c ollaboration (lAbdo et al ] l2010ch . and 
NGC 4945 ( Lenain et al.l l2010). all the other sources have not 
yet been detected at high energy. We derived upper limits on 



46 significant pulsars are reported in 'Abdo et alj ll 2010Wl . and 10 
pulsar wind nebula are discussed in Ackermann et al. ( 20111). 



their photon flux, and found that these limits are fully compati- 
ble with the expectations from the model of IPavIidou & Fieldsl 
(2001) for interactions of cosmic rays with the local interstel- 
lar medium. In the case of M33 and M83, a detection with 
Fenni/LALT should be imminent, if this model holds. The ex- 
pected fluxes for the other sources in the sample are below the 
sensitivity of ferm//LAT, even after ten years of observation. 
This is consistent with the work of Fiel ds et al. I (12010), who 
estimate that there are around five normal galaxies that are in- 
dividually detectable with Fenni/LALT. A contribution from the 
collective emission of compact objects in the host galaxies is 
found to be negligible compared to such expected emission. 

Deeper observations provided by the fermZ/LAT all sky sur- 
vey strategy will possibly allow for the detection of other nor- 
mal, non-active galaxies in the near future. 
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